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The low-energy physics of three-leg frustrated antiferromagnetic spin-S tubes in the vicinity of the 
upper critical field are studied. Utilizing the effective field theory based on the spin-wave approx- 
imation, we argue that in the intermediate-interchain-coupling regime, the ground state exhibits a 
vector chiral order or an inhomogeneous magnetization for the interchain (rung) direction and the 
low-energy excitations are described by a one-component Tomonaga-Luttinger liquid (TLL). In both 
chiral and inhomogeneous phases, the Z2 parity symmetry along the rung direction is spontaneously 
broken. It is also predicted that a two-component TLL appears and all the symmetries are restored 
in the strong-rung-coupling case. 



PACS numbers: 75.10.Jm,75.10.Pq,75.30.Kz,75.40.Cx,75.50.Ee 



i 

■4-* 



I 

S3 
O 

o 



> 
in 



o 

03 



1 

S3 
O 

o 



X 



I. INTRODUCTION 

Frustrated spin systems^ have been continuously ex- 
plored for more than five decades. Frustration is consid- 
ered as an important keyword to generate exotic, uncon- 
ventional magnetic orders, disorders and excitations in- 
cluding even spin-liquid states. Actually, frustrated sys- 
tems have provided several peculiar concepts and phe- 
nomena so far: resonating-valence-bond picture, non- 
collinear orders, symmetry-unrelated degeneracy, order- 
by-disorder mechanism, etc. 

In recent years, frustrated magnets containing four- 
spin exchanges as well as standard two-spin ones have 
been intensively studied. ? In such magnets, fascinating 
magnetic orders (nematic, chiral, dimer orders, etc.), 
which order parameter is defined by products of spin op- 
erators, are shown to be present. In a sense, these new or- 
ders are a natural consequence of the four-spin exchange 
because for such an interaction, it is possible to perform a 
mean-field approximation, S?SfS%Sf (SfSf) 



S?S 



{SfSf)(SlSf). Furthermore, it is well 
known that effects of four-spin exchanges are fairly small 
in a large number of real magnets. Thus, to discover 
intriguing magnetic orders within spin systems contain- 
ing only two-spin exchanges could more stimulate many 
experimentalists and would be theoretically a more chal- 
lenging issue. 

In one dimension, as representatives of geometrically 
frustrated spin systems with only two-spin exchanges, 
one can consider zigzag spin chains and three-leg antifer- 
romagnetic (AF) spin tubes, i.e., ladders with a periodic 
boundary condition (PBC) along the interchain (rung) 
direction. In this paper, we study the latter model in a 
magnetic field. The Hamiltonian is written as 



H = 



££[ 



JSij ■ Sij+i + J±Sij ■ Si + ij — HSfj ,(1) 



1=1 



where Sij is spin-S" operator on site j of the Zth chain 



(I = 1,2,3), J > (Jj_ > 0) is the intrachain (in- 
terchain) coupling, and the PBC S4J = Sij is im- 
posed. Focusing on the vicinity of the upper critical 
field and applying an effective field theory approach, 
we show the possibility of two interesting long-range- 
ordered states: for a certain high-magnetic-field area, a 
vector chirality (Vf •) = {(Sij x Si+ij) z ) or an inhomo- 
geneous magnetization along the rung direction occurs 
in a one-component Tomonaga-Luttinger-liquid (TLL) 
state. In the chiral phase, the Z2 rung-parity symme- 
try Sfj <-* S , ; Q ^ lj , by which Vy changes its sign, is 
spontaneously broken, while the inhomogeneous mag- 
netization in another phase breaks the one-site trans- 
lational symmetry for the rung as well as the rung- 
parity one. We also predict that a two-component TLL 
emerges and all the symmetries are preserved in the 
strong-rung-coupling regime. Recently a spin tube mate- 
rial [(CuCl2tachH) 3 Cl]Cl2 (Ref. y) has been synthesized 
and its magnetic properties could be described by a three- 
leg frustrated spin-tube modeli^^ This also promotes 
the motivation of studying the spin tube (fTj). 

Existing results of the model ([T]) are summarized here. 
In the S = \ case, the zero-field ground states are gapped 
and doubly degenerate with spontaneously breaking the 
one-site translational symmetry along the chain, at least 
when J± > 0.5J. 7 In addition, a semi-quantitative 
ground-state phase diagram in the J±-H plane (Jj_ > 0), 
which only shows gapless and gapful regimes, is con- 
structed in Ref. [H; there exists an intermediate magne- 
tization plateau with M = (Sf ■) = 1/6. In the case of 
S = integer and H = 0, the system is predicted to be 
always gapful and to conserve all symmetries.— 

Before analyzing the quantum spin tube |T]) , to discuss 
its classical version is instructive. The classical ground 
state is an umbrella structure as in Fig. [1] In this state, 
symmetries of the U(l) spin rotation around the spin 
z axis, one-site translations, and parity transformations 
along both the chain and the rung directions are all bro- 
ken. Consequently, the system exhibits a finite vector 
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spin moment 



,U(1) degeneracy 
around the S z axis 




umbrella structure 
FIG. 1: Classical ground state of the spin tube {l}. 
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From this re- 



'1,3 1 ~ 2 V 1 S' 2 (iJ+3J±y J ■ 

suit, the vector chiral order is expected to exist even in 
the quantum version. However, since generally quantum 
fluctuation is quite strong in one dimension and tends to 
destroy any ordering, it is nontrivial whether or not the 
chiral order remains and broken symmetries are restored 
in the model (fTJ) . 



II. EFFECTIVE THEORY 

Here we construct the effective theory for the quantum 
spin tube (P) in a high magnetic field. Let us begin with 
the fully polarized state with M — S. For the state, 
the energy dispersion of one magnon with AS Z = — 1 is 
exactly calculated as 

e K (k) = H-2S(J + J±) + 2SJcosk + 2SJ±cosK, (2) 

where K (= 0,±^) is the wave number for the rung 
and that for the chain, k, is in |fc| < tt. The lowest 
bands e ± 2*_ are always degenerate due to the rung-parity 
symmetry, transformation of which induces K — > —K. 
As we explain in Fig. [2j when H becomes lower than the 
upper (lower) critical value if" = AS J + 3SJ± (H l c = 
3SJ±), magnons of the lowest bands begin to condense 
(are fully condensed). Moreover, as H < H' c = AS J, 
magnons in the remaining band eo are also condensed. 

Supposing that multimagnon bound states are absent 
or their excitation energies are higher than those of one- 
magnon states (this is highly expected in antiferromag- 
netic systems and we have numerically verified it near 
the saturation), we may describe the low-energy physics 
around H ~ if" using one- magnon excitations. A suit- 
able method for such a description is spin-wave theory 
(1/5 expansion). It makes spins bosonize as 
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b} d y/2S-THj, 



(3) 



where bij is the magnon annihilation operator, and 
b\ -bi t j denotes the magnon number. Substitut- 



ing Eq. ([3]) in the model (fTJ) and introducing the Fourier 
transformation of bij for the rung as 



V K=0,±2iz/3 
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(4) 
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FIG. 2: Magnon bands in Eq. ©. 



we obtain the bosonic spin-wave Hamiltonian. As ex- 
pected, the bilinear part of bn,j reproduces the free-spin- 
wave dispersion e_R-(fc). In order to study the low-energy 
and long-distance properties of the spin tube, we further 
introduce continuous boson fields ^> q as follows: 



b 0J -> (-l) j ^Va(x), b d 



(-iy^±(x), (5) 



where an, is the lattice spacing, and x = jao. Using these 
and taking into account the magnon interaction terms up 
to the lowest order of the l/S expansion, we arrive in the 
following effective Hamiltonian, 



Ties = dx 



9=0. + ,- 



2m„ H 



PqPq 



+<Jopl + gi{p+ + P-f + foPo{p+ +P-) + fiP+P- 
+Ao(*g*5.*L +h.c.) 

+Ai(*^* ( ! ) 1' t _ + ^ 2 _^l + h.c.) + • • • , (6) 



where p q — ^q^q is the magnon-density field. (This 
Hamiltonian can also be derived via the path-integral 
approach. 10 ) The first two terms correspond to the free- 
spin- wave part, and if the chemical potential p q is posi- 
tive, the magnon ^> q is condensed^ We set /in. = ASJ—H 
and p± = p = S(AJ + 3Jj_) - H so that and 
H' c are fixed. Other parameters in Eq. © are eval- 
uated as 1/rriq = 2SJa^ (m q — to), go = 2Jao/3, 
, 9l = (4J + 3J ± )a /6, f = 8Ja /3, h = 4Ja /3, 
Ao = (8J-3J ± )a /6, and Ai = (16J-3Jj_)a /12. These 
values would be somewhat changed due to high-energy 
modes, the curvature of the dispersion, higher-order in- 
teractions, and the hard-core property of magnons ne- 
glected in the spin- wave theory. 



III. LOWEST-BAND-MAGNON CONDENSED 
STATE 

Based on the effective theory ©, we investigate the 
spin tube near saturation. In this section, we consider 
the lowest-magnon-condensed case, where /i > 0, po < 0, 
and m&x[H l c ,H' c ] < H < if". For this case, the low- 
energy physics must be governed by two condensed fields 
ty±. The effective theory is derived by integrating out 
the massive magnon "Fq. y i& the cumulant expansion in 
terms of the free-spin-wave part of 'Fo m the partition 
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function. The main effect of the V^o sector is that an 
attractive interaction between p + and p_ originates from 
the second cumulant of the Ao term. As a result, the 
coupling constant /i is changed as 

h - /i=/i-C-^=, (7) 

where C is a positive dimcnsionless constant of 0(1). 
Here, we have approximated the Matsubara Green's 
function (T T ^ (x, r)^J(0, 0)) as 1/ao (zero) when \x\ 
and JetoT are smaller (larger) than the correlation 
length (m|/io|) _1 ' 2 [t : imaginary time], and assumed 
that (m|/io|) -1 / 2 is at most O(ao).— For the resultant 
Hamiltonian TC' cS [^>±], the Haldane's harmonic- fluid ap- 
proach (i.e., bosonization) (Refs. [TMl3fT5l ) could be ap- 
plicable. Using the bosonization formulas p±(x) m 
{p ± + d x <f> ± /?r} EZ-oo e i2n ^±-*P±^ and ^ ± ~ {p ± + 

d x (j>±/ir} 1/2 Y.n=-oo el2n( ' t>± ~* P±x)e ~ iB± i where P± = 
(p±)i we obtain a bosonized Hamiltonian of the phase 

fields (</>±,#±)- Introducing the new fields 4> s ^ a = (</>+ ± 
(j)-)/V2 and 6 3>a = (0+ ±6L)/\/2 further, we can repre- 
sent the phase-field Hamiltonian as 

q—s,a 

+ 50 cos(2%/20 Q ) + 30 cos(3V2(9 Q ) + • • • , (8) 

where we have assumed p + = p- = p (see below) and 
dropped terms with spatially oscillating factors e l2nn P x . 
The g,!, and ge terms for example originate from p+p- 
and the third cumulant, respectively. Unfortunately, 
the values of g^g cannot be evaluated quantitatively 
within the present approach. In the phase-field pic- 
ture, a spin rotation around the S z axis Sj- — * e^S,-. 
the one-site translation along the chain 5*,- — > ST^+i, 
that along the rung Sj*. — > Sj+i,-) and the site-parity 
transformation along the chain — * Sf_j are, respec- 
tively, expressed as 9± — > #± + 7, (</>±(x), 0±(x)) — > 
(0±(a; + ao) - 7rp±a , 0±(x + oq) - 7r), 6>± -> 0± ± 27r/3, 
and (<£±(a;),0±(x)) — > (— 0±(— x), 0±(— a;)). Further- 
more, the rung-parity transformation Sfj «-> S^- may 
be realized by p + = ,o_ and (</>±, 9±) — > (</> T , Owing 
to these symmetries, in all vertex operators without oscil- 
lating factors, only cos[2n(0 + — (/>-)] and cos[3n(0_|_ — Q-)\ 
are allowed to exist in Eq. ([5]). The most relevant n = 1 
terms indeed appear in Eq. (jHJ). 

The bosonization approach for TL' cS evaluates the ve- 
locity v a as v a m (— fip/m) 1 / 2 . Therefore, if /1 > 0, then 
v a becomes imaginary and it means that the bosonization 
is invalid. To understand the physical meaning of this in- 
stability^ we should consider the magnon-density part 
in Ti.' cS and then define the following Ginzburg-Landau 
(GL) potential: 

T = 9i{p+ + P-f + IiP+P- -KP+ + P-)- ( 9 ) 



It is clear that as f\ > 0, the potential is minimized by 
imposing p + 7^ Moreover, it is found that 

P+ ~ P- « % d h,j - ~ b -^,;b_ ¥d ~ 2 (10) 

We thus conclude that for fx > 0, a finite long-range 
vector chiral order (Vf^) exists, and the rung-parity 
symmetry is spontaneously broken. For J± <C J (i.e., 
|/xo|/J « 1) or Jl » J [i.e., A - -0(J X )], A < 
generally holdsr^l while for J± ~ O(J) (i.e., Ao ~ 0), 
when becomes closer to ff", /1 increases and tends 
to be positive. Consequently, the chiral phase is present 
in an intermediate-rung-coupling regime. Supposing that 
p + > p_ holds in the chiral phase, we can speculate that 
the \E'_ mode constructs a massive spectrum, whereas 
the part provides a TLL stated Namely, the co- 
existence of the chiral order and the TLL is predicted. 
The presence of the TLL is also supported by the pre- 
vious study in Ref. 8. If H ~ ff", the TLL parameter 
would be close to the universal value 1. The correlation 
function of the chirality might exhibit a power decay: 
(VfjV/Jo) « (V^) 2 - const/j 2 + ... at j ^ 00^ 

Let us now discuss the case of f± < 0, where p + = p_ is 
restored and the bosonization is available. The (f> s sector 
in Eq. ([5]) yields a TLL, which is strongly stabilized by 
symmetries, while the low-energy physics of the 4> a sec- 
tor depends on whether cos(2v / 20 a ) and cos(3v / 2^a) are 
relevant or not: the scaling dimensions of these two are 
2K a and 9/(2K a ), respectively. The Hamiltonian H' cS 
leads to K a oc (—p/fi) 1 ^ 2 - Therefore, when /1 ~ and p 
is large enough, K a is always much larger than 1. At this 
case, cos(3\/20 a ) and cos(2v / 20 a ) are respectively highly 
relevant and irrelevant, and then the <p a sector obtains 
a massive spectrum. If gg > (< 0), the phase field 8 a 
is pinned on lines 9 a — y/2(2n + l)n/6 (>/2mr/3) in the 
6 + -9- plane. Among these lines, only six lines intersect 
the physically relevant "Brillouin" zone, — tt < 6+ < it 
and — 7T < < tt. This result implies that the ground 
states possess the sixfold degeneracy. To investigate the 
physical meaning of locking a and the ground-state de- 
generacy, let us focus on the magnetization per site. The 
bosonization represents it as 

(Sfj) MM- !^o( cos (V20 a + |ttZ) ) + •••. (11) 

One can see that the second term in Eq. (fTTjl causes 
a down-down-up magnetization structure in the case of 
ge > 0, while for ge < an up- up-down structure occurs: 
for instance, if 9 a is locked to zero for ge < 0, (Sfj) — 

(Sij) = M + S and (Sl d ) = M - 26 [S cx (cos(V2(9 a ))]. 
We thus conclude that an inhomogeneous magnetization 
for the rung is induced by pinning 9 a . Obviously, the 
parity and translational symmetries for the rung direc- 
tion are spontaneously broken in this state. Three of 
the sixfold degenerated states are indeed explained by 
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this inhomogeneous distribution. The meaning of the re- 
maining twofold degeneracy is unknown^ Remarkably, 
the inhomogeneously magnetized phase is not at all ex- 
pected from the classical tube system (see Fig. [1} . We 
note that this inhomogeneous distribution might slightly 
be modified if cos(3nV29 a ) with n > 2 are also rele- 
vant^ From the predictions of the chiral order for /i > 
and the inhomogeneous phase under the condition /i < 
and |/i| ~ 0, the boundary /j = is expected to be a 
first-order transition. 

When —fi/p increases so that K a < 9/4, cos(3\/2#a) 
becomes irrelevant and the low-energy physics of the 4> a 
sector is described by a Gaussian model. This transi- 
tion must be of a Beresinskii-Kosterlitz-Thouless (BKT) 
type^S After the transition, the system is in a two- 
component TLL phase with all symmetries enjoying. If 
—fi/p is further increased due to the growth of J± or 
the decrease of p, cos(2v / 20a) seems to become rele- 
vant. However, the exact results for the integrable Bose 
gas^i imply that in a one-dimensional Bose system with 
a short-range repulsive interaction, the TLL parameter 
is not usually smaller than 1 even when the interaction 
becomes extremely strong. The two-component TLL is 
hence expected to continue even when J± 3> J or p is 
small (see the Endnote [l?]). The prediction of the two- 
component TLL in the strong-rung-coupling regime is in 
agreement with a previous study applying the strong- 
rung-coupling approach to the S = § tube^ 



IV. THREE-BAND-MAGNON CONDENSED 
STATE 

Here, we consider the case where all three kinds of 
magnons v E'+ j _ i o are condensed. This situation could be 
realized under the condition of p > 0, po > 0, H l c < H < 
H' c , and J± < 4J/3. This means that the three-band- 
magnon condensed state is allowed to exist only in the 
weak-rung-coupling regime. Like Eq. let us introduce 
the GL potential for the present case as follows: 

Q = 9opI + 9i{p+ + P-f + foPo(p+ + P-) 

+fiP+P-~ MoPo- p(p++P-). ( 12 ) 

To find the stable magnon-density profile (po, p±, 
the Hessian matrix H^j — [ g®. g p ) is useful. At the lo- 
cal minimum point (poiPi P) satisfying dQ/dpj = 0, the 
eigenvalues of Hij are —4 J/3, Ci, and C2 (—4 J/3 < 
C\ < and C2 > 0). The corresponding eigenvectors are 
(8po,Sp+,5pJ) « (0,1,-1), (-C 3 ,l,l), and (C 3 ,l,l), 
where C3 > 0. The negative eigenvalue — 4J/3 and its 
eigenvector indicate that the ground state takes p+ — 
P- 7^ 0. Moreover, a positive eigenvalue C2 implies the 
existence of the TLL. We therefore predict that the chi- 
ral order (p + 7^ p_ ) and a one-component TLL state still 
remain when the system moves from the lowest-magnon- 
condensed regime to all-magnon-condensed onei 23 ' 24 At 
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FIG. 3: Schematic ground-state phase diagram of the S = 5 
spin tube (JTJ . The area away from the saturation is discussed 
elsewhere (Ref.[H). See the Endnotes and[T3. 



the boundary between these two regime, one might ob- 
serve a weak singularity such as a magnetization cusp. 



V. SUMMARY AND DISCUSSIONS 

We have studied the three-leg frustrated spin tube ([1]) 
near the upper critical field. It has been predicted that 
the vector chiral order or the inhomogeneously magne- 
tized order emerges in the magnetic-ficld-driven TLL 
phase in the intermediate-rung-coupling regime. It is 
remarkable that in these two phases, the TLL critical- 
ity (massless modes) and the spontaneous breakdown 
of discrete parity or translational symmetries for the 
rung direction coexist. We have also shown that when 
the rung coupling becomes strong enough, the inhomo- 
geneous phase vanishes and instead the two-component 
TLL occurs with preserving all the symmetries. 

Combining our results and the existent ones,— ^ we can 
draw the ground-state phase diagram for the S = 5 tube 
as in Fig. [3] The global phase structure near the sat- 
uration would common to all the cases with arbitrary 
S, as far as S < 0(1). Although in general the spin- 
wave approach used in this paper is not very reliable for 
small- S cases, we believe that it is valid if we consider 
the region where M is sufficiently close to the saturation 
value: in such a region, multimagnon scattering processes 
are expected to be negligible. When J± is changed from 
+0 to +00 with M fixed near the saturation, the fol- 
lowing scenario is expected: TLL plus chirality —►[first- 
order transition]^ TLL plus inhomogeneous magnetiza- 
tion — >[BKT transition]—* two-component TLL. 

We finally note that the predicted first-order and BKT 
transitions could not be detected by observing the mag- 
netization M because H couples to d x (j) s and p + +p_, but 
it does not directly interact 4> a and p+ — p- . A specific- 
heat measurement would be efficient in the detection. 
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